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Abstract. We present a (1 + e)-approximation algorithm running in 0(/(e) ■ n log'* n) time for finding 
the diameter of an undirected planar graph with non-negative edge lengths. 

1 Introduction 

The diameter of a graph is the largest distance between two vertices. Computing it is among the 
most fundamental algorithmic graph problems. In general weighted graphs, as well as in planar 
graphs, the only known way to compute the diameter is to essentially solve the (more general) All- 
Pairs Shortest Paths (APSP) problem and then take the pair of vertices with the largest distance. 

In general weighted graphs, solving APSP (and thus diameter) currently requires O(n^) time. 
The fastest algorithm to date is 0(n'^(loglogn)^/log^n) by Chan [5], or for sparse graphs 0{mn + 
v? logn) by Johnson [12], with a small improvement to 0{mn + n? log log n) in [18]. 

In weighted p/anar graphs, solving APSP can be done in 0{'n?) time by Frederickson [10]. While 
this is optimal for APSP, it is not clear that it is optimal for diameter. Currently, only a logarithmic 
factor improvement by Wulf-Nilsen [20] is known for the diameter, running in 0(n^ (log log n)^/ log n) 
time. A long standing open problem [6] is to find the diameter in truly subquadratic 0{n'^^^) time. 
Eppstein [8] has shown that if the diameter in a planar graph is bounded by a fixed constant then 
it can be found in 0{n) time. Fast algorithms are also known for some simpler classes of graphs 
like outer- planar graphs [9], interval graphs [17], and distance-hereditary graphs [7]. 

In lack of truly subcubic-time algorithms for general graphs and truly subquadratic-time algo- 
rithms for planar graphs it is natural to seek faster algorithms that approximate the diameter. It is 
easy to approximate the diameter within a factor of 2 by simply computing a Single-Source Shortest 
Path (SSSP) tree from any vertex in the graph and returning twice the depth of the deepest node in 
the tree. This requires 0{m + n\ogn) time for general graphs and 0(n) time for planar graphs [11]. 
For general graphs, Aingworth et al. [2] improved the approximation factor from 2 to 3/2 at the 
cost of 0{m^/n + 'n?) running time, and Boitmanis et al. [4] gave an additive approximation factor 
of 0(\/n) with 0{my/n) running time. For planar graphs, the current best approximation is a 
3/2-approximation by Berman and Kasiviswanathan running in 0(^3/2) time [3]. We improve this 
to a (1 -|- e)-approximation running in 0{n) time for any fixed < e < 1. 

Our Result. 

Theorem 1. Given an undirected planar graph with non-negative edge lengths and diameter d, 
for any e > we can compute an approximate diameter d' (where d < d' < (1 + e) ■ d) in time 
0(n log^ n/e4 + n-20(i/")). 



Summary of the Algorithm. A lemma of Lipton and Tarjan [15] states that, for any SSSP tree 
T in a planar graph, there is a non-tree edge e (where e might possibly be a non-edge of the planar 
graph) such that the strict interior and strict exterior of the unique simple cycle C in TU {e} each 
contains at most 2/3 • n vertices. The vertices of C therefore form a separator consisting of two 
shortest paths with the same common starting vertex. 

Let Gin (resp. Gout) be the subgraph of G induced by G and all interior (resp. exterior) vertices 
to G. Let d{Gin,Gout,G) denote the largest distance in the graph G between a marked vertex in 
ViGin) and a marked vertex in V{Gout)- In the beginning, all vertices of G are marked and we seek 
the diameter which is d{G,G,G). We use a divide and conquer algorithm that first approximates 
d{Gin,Gout,G), then unmarks all vertices of C, and then recursively approximates d{Gin,Gin,G) 
and d{Gout, Gout, G) and takes the maximum of all three. We outline this algorithm below. Before 
running it, we compute an SSSP tree from any vertex using the linear-time SSSP algorithm of 
Henzinger et al. [11]. The depth of the deepest node in this tree already gives a 2-approximation 
to the diameter d{G, G, G). Let x be the obtained value such that x < d{G, G, G) < 2x. 

Reduce d{Gin, Gout, G) to d{Gin, Gout, Gt) in a tripartite graph Gt-' The separator G is composed of 
two shortest paths P and Q emanating from the same vertex, but that are otherwise disjoint. We 
carefully choose a subset of 16/e vertices from G called portals. The first (resp. last) 8/e portals are 
evenly spread vertices across the prefix of P (resp. Q) of length 8x. The purpose of the portals is to 
approximate a shortest u-to-v path for u G Gin and v G Gout by forcing it to go through a portal. 
Formally, we construct a tripartite graph Gt with vertices (V {G in), portal s,V (Gout))- The length 
of edge {u G V{Gin),v G portals) or {u £ portals, v £ V{Gout)) in Gt is the u-to-v distance in G. 
This distance is computed by running the SSSP algorithm of [11] from each of the 16/e portals. By 
the choice of portals, we show that d{Gin, Gout, Gt) is a (1 + 2e)-approximation of d{Gin, Gout, G). 

Approximate d{Gin, Gout, Gt): Hi is the maximum edge-length of Gt, then note that d{Gin, Gout, Gt) 
is between d. and 2L This fact makes it possible to round the edge- lengths of Gj to be in {1, 2, . . . , 1/e} 
so that ei ■ d{Gin,Gout,Gt) after rounding is a (1 + 2e)-approximation to d{Gin,Gout,Gt) before 
rounding. Using the fact that after rounding d(Gj„, Gout, Gt) is constant, we give a linear-time algo- 
rithm to compute it exactly, thus approximating d{Gin, Gout, G). We then unmark all vertices of C 
and move on to recursively approximate d{Gin, Gin, G) (the case of d{Gout, Gout, G) is symmetric). 

Reduce d{Gin,Gin,G) to d{Gin,Gin,Gf^) in a planar graph Gf^ of size at most 2/3 • n; In or- 
der to apply recursion, we construct a planar graph Gf^ that has at most 2/3 • n vertices and 
d{Gin,Gin,Gf^) is a (1 + e/(21ogn))-approximation^ to d{Gin,Gin,G). To construct Gf^, we first 
carefully choose a subset of 2561ogn/e vertices from G called dense portals. We then compute 
all 0((2561ogn/e)^) shortest paths in Gout between dense portals. The graph B' obtained by the 
union of all these paths has at most 0((256 log n/e)^) vertices of degree > 2. We contract vertices 
of degree = 2 so that the number of vertices in B' decreases to 0((2561ogn/e)^). Appending this 
small graph B' (after unmarking all of its vertices) as an exterior to Gin results in a graph Gf^ that 
has \Gin\ + 0((2561ogn/e)'^) vertices and d{Gin,Gin,Gf^) is a (1 + e/(21ogn))-approximation of 

in, Gin, G). 

The problem is still that the size of Gf^ is not necessarily bounded by 2/3 • n. This is because C 
(that is part of Gf^) can be as large as n. We show how to shrink to size roughly 2/3 • n while 

^ log n = logj n throughout the paper. 
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d{Gin,Gin,Gf^) remains a (1 + e/(2 log n))-approximation of d{Gin,Gin,G). To achieve this, we 
shrink the G part of Gf^ so that it only includes the dense portals without changing d{Gin, Gin, Gf^). 

Approximate d{Gin, Gin, Gf^): Finally, once < 2/3 • n we apply recursion to d{Gin, Gin, ^m)- 
In the halting condition, when |G^| < (2561ogn/e)'^, we naively compute d{Gin,Gin,Gf^) using 
APSP. 

Related Work. The use of shortest-path separators and portals to approximate distances in 
planar graphs was first suggested in the context of approximate distance oracles. These are data 
structures that upon query u,v return a (1 + e)-approximation of the u-to-v distance. Thorup [19] 
presented an 0(l/e • nlogn)-space oracle answering queries in 0(l/e) time on directed weighted 
planar graphs. Independently, Klein [14] achieved these same bounds for undirected graphs. 

In distance oracles, we need distances between every pair of vertices and each vertex is associated 
with a possibly different set of portals. In our diameter case however, since we know the diameter is 
between x and 2x, it is possible to associate all vertices with the exact same set of portals. This fact 
is crucial in our algorithm, both for its running time and for its use of rounding. Another important 
distinction between our algorithm and distance oracles is that distance oracles upon query {u, v) 
can inspect all recursive subgraphs that include both u and v. We on the other hand must have 
that, for every (u, v), the shortest u-to-v path exists (approximately) in the unique subgraph where 
u and V are separated by C. This fact necessitated our construction of Gf^ and G^^^. 

2 The Algorithm 

In this section we give a detailed description of an algorithm that approximates the diameter of an 
undirected weighted planar graph Q = (V, E) in the bounds of Theorem 1. The algorithm computes 
a (1 + e)-approximation of the diameter d = d{G, G, G) for G = Q. This means it returns a value 
d' where d < d' < {1 + e) ■ d (recall that, before running the algorithm, we compute a value x such 
that X < d < 2x hy computing a single-source shortest-path tree from an arbitrary vertex in ^). We 
focus on approximating the value of the diameter. An actual path of length d' can be found in the 
same time bounds. For simplicity we will assume that shortest paths are unique. This can always 
be achieved by adding random innitesimal weights to each edge; using the isolation lemma of [16]. 
Also, to simplify the presentation, we assume that e < 0.1 and we describe a (l + 7e)-approximation. 
(then just take e' = e/7). 

The algorithm is recursive and actually solves the more general problem of finding the largest 
distance only between all pairs of marked vertices. In the beginning, we mark all n = \V{G)\ 
vertices of G = G and approximate d{G, G, G) (the largest distance in G between marked vertices 
in V{G)). Throughout the recursive calls, we maintain the invariant that the distance between any 
two marked vertices in the graph G of the recursive call is a (1 + e)-approximation of their distance 
in the original graph Q (there is no guarantee on the marked-to-unmarked or the unmarked-to- 
unmarked distances). We denote by 6g{u, v) the u-to-v distance in the original graph Q. 

The recursion is applied according to a variant of the shortest-path separator decomposition 
for planar graphs by Lipton and Tarjan [15]: We first pick any marked vertex t^i and compute in 
linear time the SSSP tree from vi in G. In this tree, we can find in linear time two shortest paths 
P and Q (both emanating from vi ) such that removing the vertices oi G = PUQ from G results in 
two disjoint planar subgraphs A and B (i.e., there are no edges in V{A) x V{B)). The number of 
vertices of G can be as large as n but it is guaranteed that |T^(j4)| < 2/3 • n and |V^(S)| < 2/3 • n. 
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Notice that the paths P and Q might share a common prefix. It is common to not include this 
shared prefix in C. However, in our case, we must have the property that P and Q start at a marked 
vertex. So we include in C the shared prefix as well (See Fig. 1). 




Fig. 1. Shortest-path separator: In the above weighted undirected planar graph G, the black nodes 
constitute the shortest path separator C composed of two shortest paths P and Q emanating from 
the same vertex v\ . The subgraph of G induced by the white nodes is denoted A and the subgraph 
of G induced by the gray nodes is denoted B. The graph Gin is the subgraph induced by all white 
and black vertices and the graph Gout is the subgraph induced by all gray and black vertices. 

Let Gin (resp. Gout) be the subgraph of G induced by V{C) U V{A) (resp. V{G) U V{B)). 
To approximate d{G,G,G), we first compute a (1 + 5e)-approximation di of d{Gin,Gout,G) (the 
largest distance in G between the marked vertices of V{Gin) and the marked vertices of V{Gout))- 
In particular, di takes into account all V{G) x V{G) distances. We can therefore unmark all the 
vertices of C and move on to approximate d2 = d{Gin, Gin, G) (approximating ^3 = d{Gout, Gout, G) 
is done similarly). We approximate d{Gin, Gin, G) by applying recursion on d{Gin, Gin, Gfu) where 
|y(G^)| < 2/3 • n. The marked vertices in Gf^ and in Gin are the same and d{Gin,Gin,Gf^) is a 
(1 + e/(21ogn))-approximation of d{Gin,Gin,G). This way, we add a factor of e/(21ogn) to the 
diameter in each recursive call. Since the recursive depth is O(logn) (actually, it is never more than 
1.8 logn) we get a (1 + 5e) • (1 + e) < (1 + 7e)-approximation d2 to d{Gin, Gin, G). Finally, we return 
the maximum of di, (i2, ^3. 

2.1 Reduce d{Gin, Gout, G) to d{Gin, Gout, Gt) 

Our goal is now to approximate d{Gin,Gout,G). For u G Gin and v G Gout, we approximate a 
shortest u-io-v path in G by forcing it to go through a portal. In other words, consider a shortest 
u-io-v path. It is is obviously composed of a shortest u-to-c path in G concatenated with a shortest 
c-io-v path in G for some vertex c G C. We approximate the shortest u-to-v path by insisting that 
c is a portal. The fact that we only need to consider u-io-v paths that are of length between x and 
2x makes it possible to choose the same portals for all vertices. 

We now describe how to choose the portals in linear time. Recall that the separator G is 
composed of two shortest paths P and Q emanating from the same marked vertex vi . The vertex 
vi is chosen as the first portal. Then, for i = 2, . . . ,8/e we start from Vi^i and walk on P until we 
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reach the first vertex whose distance from Vi-i via P is greater than ex. We set this vertex as the 
portal Vi and continue to i + Notice that P is a shortest path in G and the portals are chosen in 
a prefix of P of length 8x. This might seem counterintuitive as we know that any shortest path P 
in the original graph Q is of length at most 2x. However, since one endpoint of P is not necessarily 
marked, it is possible that P is a shortest path in G but not even an approximate shortest path in 
the original graph Q. We do the same for Q, and we get a total of 16/e portals (See Fig. 2). 



Fig. 2. Portals: The six circled vertices are the 16/e portals in the 8x prefixes of P and Q. The 
shortest path between u and v goes through the separator vertex c and is approximated by the 
u-to-p{c) and the p{c)-to-v shortest paths where p{c) is the closest portal to c. The distance from 
c to p(c) is at most ex. 

Once we have chosen the portals, we move on to construct a tripartite graph Gt whose three 
vertex sets are {V{Gin),portals,V{Gout))- The length of edge (n G V{Gin),v G portals) or (n S 
portals, V G V{Gout)) is the u-to-v distance in G. This distance is computed by running the linear- 
time SSSP algorithm of Henzinger et al. [11] in G from each of the 16/e portals in total 0{l/e ■ 
\V{G)\) time. The following lemma states that our choice of portals implies that d{Gin, Gout, Gt) is 
a good approximation of d{Gin, Gout, G). 



Lemma 1. Ifd{Gin, Gout, G) > x, then d{Gin, Gout, Gt) is a {l+2e) -approximation ofd{Gin, Gout, G) 
If d{Gin, Gout, G) < X then d{Gin, Gout, Gt) < (1 + 2e) ■ x. 



Proof. The first thing to notice is that d{Gin,Gout,Gt) > d{Gin,Gout,G). This is because every 
shortest u-to-v path in Gt between a marked vertex u E V{Gin) of the first column and a marked 
vertex v G V{Gout) of the third column corresponds to an actual u-to-v path in G. 

We now show that d{Gin, Gout, Gt) < (1 + 2e) • d{Gin, Gout, G). We begin with some notations. 
Let Pt denote the shortest path in Gt realizing d{Gin,Gout,Gt). The path Pt is a shortest u-to-v 
path for some marked vertices u E Gin and v S Gout- The length of the path Pt is 6ctiu,v). Let 
Pq denote the shortest u-to-v path in G that is of length Sciujv) and let Pg denote the shortest 
u-to-v path in the original graph Q that is of length 5g{u, v). Recall that we have the invariant that 
in every recursive level for every pair of marked vertices 5g{'^, v) < {1 + e) ■ 5g{u, v). We also have 
that 6g{u, v) < 2x and so 6g{u, v) < 2x ■ {1 + e). For the same reason, since vi (the first vertex of 
both P and Q) is also marked, we know that 5g{vi,u) is of length at most 2x ■ {1 + e). 
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The path must include at least one vertex c £ C. Assume without loss of generality that 
c E P. We claim that c must be a vertex in the prefix of P of length 8x. Assume the converse, 
then the vi-to-c prefix of P is of length at least 8x. Since P is a shortest path in G, this means 
that 6c{vi,c) is at least 8x. However, consider the fi-to-c path composed of the vi-to-u shortest 
path (of length 6c{vi,u) < 2x • (1 + e)) concatenated with the u-to-c shortest path (of length 
6g{u,c) < Sciujv) < 2x • (1 + e)). Their total length is 4x • (1 + e) which is less than 8x (since 
e < 1) thus contradicting our assumption. 

After establishing that c is somewhere in the 8x prefix of P, we now want to show that (u, v) < 
(1 + 2e) • 6g{u,v). Let p(c) denote the closest portal to c on the path P. Notice that by our 
choice of portals and since c is in the 8x prefix of P we have that 5g{c,p{c)) < ex. By the 
triangle inequality we know that (5g(u,p(c)) < 5g{u,c) + 6g{c,p{c)) < 5Giu,c) + ex and similarly 
^g{p{c)iv) < 5g{c,v) + ex. This means that 

d{Gin,Gout,Gt) = 6Gt{u,v) 

< 6g{u,p{c)) + 5g(p(c),w) 

< 6g{u, c) + (5g(c, v) + 2ex 
= 6g{u, v) + 2ex 

< d{Gin,Gout,G) + 2ex 

< il + 2e)-d{Gin,GouuG). 

In the last inequality we assumed that d{Gin, Gout, G) > x. If d{Gin, Gout, G) < x, then we get that 
d{Gin, Gout, Gt) < (1 + 2e) ■ x. The lemma follows. □ 

By Lemma 1, approximating d{Gin, Gout, G) when d{Gin, Gout, G) > x reduces to approximating 
d{Gin,Gout,Gt). The case of d{Gin,Gout,G) < x means that the diameter d of the original graph 
G is not a (u £ Gin)-to-{v G Gout) path. This is because d > x > d{Gin, Gout, G) > d{Gin, Gout,Q)- 
So d will be approximated in a different recursive call (when the separator separates the endpoints 
of the diameter). In the meanwhile, we will get that d{Gin, Gout, Gt) is at most (1 + 2e) ■ x and so 
it will not compete with the correct recursive call when taking the maximum. 

2.2 Approximate <i(Gin, Gout, Gt) 

In this subsection, we show how to approximate the diameter in the tripartite graph Gt- We give 
a (1 + 2e)-approximation for d{Gin,Gout,Gt). By the previous subsection, this means we have a 
(1 + 2e)(l + 2e) < (1 + 5e)-approximation for d{Gin, Gout, G). From the invariant that distances in 
G between marked vertices are a (1 + e)-approximation of these distances in the original graph 
we get a (1 + 5e)(l + e) < (1 + 7e)-approximation for d{Gin, Gout, Q) in the original graph Q. 

We now present our (1 + 2e)-approximation for d{Gin,Gout,Gt) in the tripartite graph Gt- 
Recall that Pt denotes the shortest path in Gt realizing d{Gin,Gout,Gt). By the definition of Gt, 
we know that the path Pt is composed of only two edges: (1) edge {u,p) between a marked vertex 
u of the first column (i.e., u £ V{Gin)) and a vertex p of the second column (i.e., p corresponds 
to some portal in G). (2) edge {p,v) between p and a marked vertex v of the third column (i.e., 
V £ ViGout))- 

Let X (resp. Y) denote the set of all edges in Gt adjacent to marked vertices of the first 
(resp. third) column. Let i denote the maximum edge-length over all edges in X L)Y. Notice that 
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I < d{Gin, Gout, Gt) < 2^. We round up the lengths of ah edges in X U y to the closest multiple 
of el. The rounded edge- lengths are thus all in {e^, 2e^, 3e£, ...,£}. We denote Gt after rounding 
as G[. Notice that d{Gin, Gout, G'l) is a (1 + 2e)-approximation of d{Gin, Gout, Gt)- This is because 
the path Pt is of length at least i and is composed of two edges, each one of them has increased its 
length by at most el. 

We now show how to compute d{Gin,Gout,G[) exactly in linear time. We first divide all the 
edge-lengths of Gj by el and get that G'l has edge-lengths in {1, 2, 3, . . . , 1/e}. After finding 
d{Gin,Gout,G'i) (which is now a constant) we simply multiply the result by el. The following 
lemma states that when the diameter is constant it is possible to compute it exactly in linear time. 
Note that we can't just use Eppstein's [8] linear-time diameter algorithm for a planar graph whose 
diameter is bounded by a fixed constant since in our case we get a non-p/anar tripartite graph G[. 

Lemma 2. d{Gin,Gout,G[) can be computed exactly in time 0{\V{G)\/e + 2^'^^^^^). 

Proof. Recall that in G'^ we denote the set of all edges adjacent to marked vertices of the first and 
third column as X and Y . The length of each edge in X U y is in {1, 2, . . . , The number of 
edges in X (and similarly in Y) is at most 16|y(G)|/e. This is because the first column contains 
\Gin\ < |^(G)| vertices and the second column contains 16/e vertices vi,V2, ■ ■ ■ , ^le/e (the portals). 

For every marked vertex v in the first (reps, third) column, we store a 16/e-tuple vx (reps, vy) 
containing the edge lengths from v to all vertices of the second column. In other words, the tuple 
Vx = vi), 5{v, V2), . . ■ , S{v, fi6/e)) where 6{v, Vi) is the length of the edge {v, Vi). 

Notice that every 5{v,Vi) is in {1, 2, . . . , Furthermore, every |(5(f,z;j+i) — 5{v.,Vi)\ is in 

{0, 1, 2}. To see this, notice that vi and fj+i are adjacent portals on the separator (P or Q). So the 
distance between them in G is 5{vi,ViJ^i) < ex, and since we scaled (divided by ei) 5{yi, viJ^i) is at 
most ex /el. Because x < 21 we get that ex/el < 2. Finally, since 5{v,Vi) and 6{v,Vi-\-i) correspond 
to distances, by the triangle inequality we have \5{v,Vi+i) — 5{v,Vi)\ < 6{vi,Vi+i) < 2. Overall, we 
get that a tuple vx or vy has 16/e entries, the first entry is 6{v, vi) £ {1,2,..., 1/e} and for every 
other entry i + 1 it holds that 6{v, t'i+i) — S{v, Vi) can be one of 5 values: —2, —1, 0, 1, 2. The total 
number of tuples is 16|y(G)|/e but the number of different tuples is therefore only k = (1/e) •5^^/''. 

We create two binary vectors Vx and Vy each of length k. The i'th bit of Vx (reps. Vy) is 1 
iff the z'th possible tuple exists as some vx (reps. vy). Creating these vectors takes 0{\V{G)\/e) 
time. Then, for every 1 bit in Vx (corresponding to a tuple of vertex u in the first column) and 
every 1 bit in Vy (corresponding to a tuple of vertex v in the third column) we compute the u-to-v 
distance in G[ using the two tuples in time 16/e. We then return the maximum of all such {u,v) 
pairs. Notice that a 1 bit can correspond to several vertices that have the exact same tuple. We 
arbitrarily choose any one of these. There are k entries in Vx and k entries in Vy so there are 0{k'^) 
pairs of 1 bits. Each pair is examined in 0(16/e) time for a total of 0((l/e) • k'^) = 2*^^^/^) time. □ 

To conclude, we have so far seen how to obtain a (1 + 5e)-approximation for d{Gin,Gout,G) im- 
plying a (1 + 7e)-approximation for d{Gin,Gout,Q) in the original graph Q. The next step is to 
unmark all vertices of G and move on to recursively approximate d{Gin,Gin,G) (approximating 
d{G out. Gout, G) is done similarly). 

2.3 Reduce d{Gir^, Girr, G) to Gir^, Gfj 

In this subsection we show how to recursively obtain a (1 + 5e)-approximation of d{Gin,Gin,G) 
and recall that this implies a (1 + 7e)-approximation of d(Gin, Gin, G) in the original graph Q since 
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we will make sure to maintain our invariant that, at any point of the recursion, distances between 
marked vertices are a (1 + e)-approximation of these distances in the original graph Q. 

It is important to note that our desired construction can be obtained with similar guarantees 
using the construction of Thorup [19] for distance oracles. However, we present here a simpler 
construction than [19] since, as apposed to distance oracles that require all-pairs distances, we can 
afford to only consider distances that are between x and 2x. 

There are two problems with applying recursion in order to solve d{Gin,Gin,G). The first is 
that can be as large as and we need it to be at most 2/3 • We do know 

however that the number of marfced vertices in V{Gin) is at most 2/3- |y(G)| . The second problem 
is that it is possible that the u-to-v shortest path in G for u,v £ Gin includes vertices of Gout- 
This only happens if the u-to-v shortest path in G is composed of a shortest u-to-p path {p S P) 
in Gin, a shortest p-to-q path {q £ Q) in Gmit-, and a shortest q-to-v path in Gin- To overcome 
these two problems, we construct a planar graph Gf^ that has at most 2/3 • |T^(G')| vertices and 
d{Gin,Gin,GfJ is a (1 + e/(21ogn))-approximation to d{Gin,Gin,G). 

Recall that the subgraph B of G induced by all vertices in the strict exterior of the separator G 
is such that < 2/3- |V^(G)| and Gout = BUC. The construction of Gf^ is done in two phases. In 
the first phase, we replace the B part of G with a graph B' of polylogarithmic size. In the second 
phase, we contract the C part of G to polylogarithmic size. 

Phase I: replacing B with B' . To construct Gf^, we first choose a subset of 2561ogn/e vertices 
from C called dense portals. The dense portals are chosen similarly to the regular portals but there 
are more of them. The marked vertex vi (the first vertex of both P and Q) is chosen as the first 
dense portal. Then, for i = 2, . . . , 1281ogn/e we start from and walk on P until we reach the 
first vertex whose distance from Vi-i via P is greater than ex/(161ogn). We set this vertex as the 
dense portal Vi and continue to i + We do the same for Q, for a total of 256 log n/e dense portals. 

After choosing the dense portals, we compute all 0((2561ogn/e)^) shortest paths in Gout be- 
tween dense portals. This can be done using SSSP from each portal in total 0{\V{Gout)\ • logn/e) 
time. It can also be done using the Multiple Source Shortest Paths (MSSP) algorithm of Klein [13] 
in total 0{\V{Gout)\ • logn + log^ n/e^) time. 

Let B' denote the graph obtained by the union of all these dense portal to dense portal paths in 
Gout- Notice that since these are shortest paths, and since we assumed shortest paths are unique, 
then every two paths can share at most one consecutive subpath. The endpoints of this subpath 
are of degree > 2. There are only 0((2561ogn/e)^) paths so this implies that the graph B' has at 
most 0((2561ogn/e)^) vertices of degree > 2. We can therefore contract vertices of degree = 2. 
The number of vertices of B' then decreases to 0((256 log n/e)^), it remains a planar graph, and 
its edge lengths correspond to subpath lengths. 

We then unmark all vertices of B' and append B' to the infinite face of Gin- In other words, we 
take the disjoint union of Gin and B' and identify the dense portals of Gin with the dense portals 
of B' . This results in a graph Gf^ that has + 0((2561ogn/e)^) vertices. In Lemma 3 we 

will show that d{Gin, Gin, Gf^) can serve as a ( 1 + ^/(2 log n))-approximation to d{Gin, Gin, G). But 
first we will shrink so that the number of its vertices is bounded by 2/3 • 

Phase II: shrinking The problem with the current Gf^ is still that the size of V{Gf^) is 

not necessarily bounded by 2/3 • |y(G)|. This is because G (that is part of V{Gf^)) can be as large 
as n. We now show how to shrink y(G^) to size 2/3 • |T^(G)| while d{Gin,Gin,Gf^) remains a 
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Fig. 3. On the left: The graph before shrinking. The white vertices are the vertices of A, the 
black vertices are the vertices of C that are not dense portals, the six red circled vertices are the 
2561ogn/e dense portals, and the gray vertices are the vertices B' \C with degree > 2. On the 
right: The graph after shrinking. The edges adjacent to vertices of C that are not dense portals 
are now replaced with edges to dense portals. 

(1 + e/(2 log n))-approximation of d{Gin,Gin,G). To achieve this, we shrink the C part of V{Gf^) 
so that it only includes the dense portals. We show how to shrink P, shrinking Q is done similarly. 

Consider two dense portals Vi and Uj+i on P (i.e., Vi is the closest portal to Vi+i on the path 
P towards vi). We want to eliminate all vertices of P between Vi and Wj+i. Denote these vertices 
by pi, . . . ,pk- If Vi is the last portal of P (i.e., i = 1281ogn) then pi, . . . ,pk are all the vertices 
between Vi and the end of P. Recall that A is the subgraph of G induced by all vertices in the strict 
interior of the separator G. Fix a planar embedding of Gf^. We perform the following process as 
long as there is some vertex u in QU A which is a neighbor of some pj , and which is on some face 
of the embedding that also contains Vi. We want to "force" any shortest path that goes through 
an edge {u,pj) to also go through the dense portal Vi. To this end, we delete all such edges {u,pj), 
and instead insert a single edge {u,Vi) of length mmj{i{u,pj) + 6c{pj,Vi)}. Here, i{u,pj) denotes 
the length of the edge {u,pj) (it may be that £{u,pj) = oo if {u,pj) is not an edge) and 6G{pj,Vi) 
denotes the length of the pj-to-Vi subpath of P. It is important to observe that the new edge {u, Vi) 
can be embedded while maintaining the planarity since we have chosen u to be on the same face 
as Vi. Observe that once the process ends, the vertices pj have no neighbors m. Q[J A. 

Finally, we replace the entire Uj+i-to-Uj subpath of P with a single edge (fi+i, Vi) whose length 
is equal to the entire subpath length. If Vi is the last dense portal in P then we simply delete the 
entire subpath between Vi and the end of P. The entire shrinking process takes only linear time in 
the size of |y(G)| since it is linear in the number of edges of G^^ (which is a planar graph). 

The following Lemma asserts that after the shrinking phase d{Gin,Gin,Gf^) can serve as a 
(1 + e/(21ogn))-approximation to d{Gin,Gin,G). 

Lemma 3. d{Gin, Gin, G) < d{Gin, Gin, ) < d{Gin, Gin, G) + ex/(21ogn) 

Proof. First observe that d{Gin,Gin,Gf^) > d{Gin,Gin,G). This is because every vertex of Gj„ 
that is marked in G is also a marked vertex in G^, and any shortest u-to-v path in Gf^ corresponds 
to an actual u-to-v path in G. 
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We now show that d{Gin, Gin, Gf^) < d{Gin, Gin, G) + ex/{2 logn). Let P+ denote the shortest 
u-to-v path in G^ reahzing d{Gin,Gin,G^^^). Both u and v are marked vertices in Gin and the 
length of is 6p+ {u,v). Let Pq denote the shortest u-to-v path in G that is of length dG{u,v). 

in 

Case 1: If Pq does not include any vertex of C then Pq is also present in and therefore 

d{Gin, Gin, < d[Gin, Gin, G). 

Case 2: If Pq includes vertices that are not in Gin (i-e., vertices in Gout \ G) then Pq must be 
composed of a shortest u-io-p path [p G P) in Gin, a shortest p-io-q path (g G Q) in Gout, and a 
shortest g-to-'u path in Gin- 

We first claim that p must be a vertex in the prefix of P of length 8x (a similar argument holds 
for q and Q). Assume the converse, then the prefix of P from vi (the first vertex of both P and Q) 
to p is of length at least 8x. Recall that wc have the invariant that in every recursive level for every 
pair of marked vertices 5g{u, v) < (1 + e) ■ dg{u, v) < 2x- (l + e). For the same reason we know that 
Sg{vi, u) < 2.x • (l+e). Since P is a shortest path in G, this means that 5g{vi,p) > 8x. However, 
consider the vi-to-p path composed of the vi-to-u shortest path (of length 5g{vi,u) < 2x • (1 + e)) 
concatenated with the u-io-p shortest path (of length 6g{u,p) < 5g{u,v) < 2x • (1 + e)). Their total 
length is 4x • (1 + e) which is less than Sx (since £ < 1) thus contradicting our assumption. 

We now show that 5(~,+ {u,v) < Sg{u,v) + ex/ (2 logn). For c G P (resp. c e Q), let p{c) denote 

in 

the first dense portal encountered while walking from c towards vi on the path P (resp. Q). Notice 
that since p and q are in the 8x prefixes of P and Q we have that SG{p,pip)) < ex/(161ogn) and 
Sg{(1,p{q)) < ex/(161ogn). From the shrinking phase, it is easy to see that G^ includes a ■u-to- 
p{p) path of length Sg{u,p) + Sg{p,p{p)) and so Sq+ {u,p{p)) < Sg{u,p) + ex/(161ogn). Similarly, 

in 

(^(9)5'") ^ Sg{q,v) + ex/(161ogn). Furthermore, since Gf was appended with shortest paths 

in 

between dense portals in Gout wc have {p{p),p{q)) < 5Gout{p{p)'p) + ^Gout{P^Q)+^GouM^Pi(l)) = 

in 

5g{p{p),p) + ^GoutiP^l) + ^G{q,p{q)) < 6G„^tip,q) + ex/ {Slog n). To conclude we get that 
d{Gin, Gin, GfJ = 5q+ {u, v) 

< 5q+ {u,p{p)) + 5q+ {p{p),p{q)) + 5q+ {p{q),v) 

%7i in %n 

< 5g{u,p) + 5G,^t{p,q) + 5G{q,v) + ex/ {Alogn) 
= 6g{u, v) -\-ex/ (4 log n) 

< d{Gin, Gin, G) + ex/(4 log n) 

< d{Gin, Gin, G) + ex/ {2 log n). 

Case 3: Finally, we need to consider the case where Pq includes only vertices of Gin- We assume 
Pg includes vertices of P and\or vertices of Q (otherwise this was handled in Case 1). We focus on 
the case that Pg includes vertices of both P and Q. The case that Pg includes vertices of one of P 
or Q follows immediately using a similar argument. 

Since P and Q are shortest paths, then Pg must be composed of the following shortest paths: 
a u-to-p path {p G P) in Gin, a p-io-p' subpath [p' G P) of P, a p'-to-q' path (q' G Q) in Gin, a 
g'-to-g subpath {q G Q) of Q, and a q-to-v path in Gin- Following the same argument as in Case 2, 
we know that p and p' (reps, q and q') must in the prefix of P (reps. Q) of length 8x. This means 
Sg{c,p{c)) < ex/(161ogn) for every c G {p,p',q,q'}- 

From the shrinking phase, it is easy to see that includes a u-to-p{p) path of length 
5g{u,p) + 6Gip,p{p)) and so 6(.+ {u,pip)) < Sg{u,p) + ex/(161ogn). Similarly, 6r,+ ip{q),v) < 
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6G{q,v)+ex/{mogn), and 5^+ {p{p'),p{q')) < SG{p{p'),p') + 5G{p',q')+SG{q',p{q')) < SG{p',q') + 

in 

ex/ (8 log n). Furthermore, since subpaths of P in Gf^ between dense portals capture their exact 
distance in G we have that 6^+ {p{p),p{p')) < 5g{p{p)jP) + ^g{p,p') + ^g{p',p{p')) < ^g{p,p') + 

in 

6x/(81ogn) and similarly S^-\- {p{q')->p{q)) < Sciq^Q) + ex/{8logn). To conclude we get that 

in 

d{Gin, Gin, Gfj = 5^+ {u, v) 

in 

< 5fj+ {u,p{p)) + 5(.+ {p{p),p{p')) + 6q+ {p{p),p{q)) + {p{q),p{q)) + 5^+ {p{q),v) 

in in in in in 

< Sg{u,p) +6g{p,p) + 6G{p,q) +6G{q,q) + 6G{q,v) +ex/(21ogn) 
= Sg{u,v) +ex/(21ogn) 

< d{Gin, Gin, G) + ex/ {2logn). 

□ 

Corollary 1. If d{Gin,Gin,G) > x then d{Gin,Gin,Gf^) is a {I + s/{2 log n))- approximation of 
d{Gin,Gin,G). If d{Gin,Gin,G) < X, then we get that d{Gin,Gin,Gfj < (1 + e/(2 log n)) • x. 

By the above corollary, approximating d{Gin, Gin, G) when d{Gin, Gin, G) > x reduces to approxi- 
mating d{Gin, Gin, GfjJ. When d{Gin, Gin, G) < x ii means that the diameter of the original graph 
Q is not a (u G Gj„)-to-(u G Gin) path and will thus be approximated in a different recursive call. 

Finally, notice that indeed we maintain the invariant that the distance between any two marked 
vertices in the recursive call to Gf^ is a (1 + e)-approximation of the distance in the original graph 
Q. This is because, by the above corollary, every recursive call adds a 1 + e/(21ogn) factor to the 
approximation. Each recursive call decreases the input size by a factor of (2/3 + o(l))^^. Hence, 
the overall depth of the recursion is at most \ogi^_^(i-^ n < 1.8 log n. Since 

(l + e/(21ogn))i-8i°s" < e^-^' < 1 + e 

the invariant follows (we assume in the last inequality that e < 0.1). Together with the (1 + 5e)- 
approximation for d{Gin,Gout,Q) in the original graph Q, we get a (1 + 5e) • (1 + e) < (1 + 7e)- 
approximation of d{Gin, Gin, G) in the original graph Q, once we apply recursion to d{Gin, Gin, Gf^)- 
We note that our recursion halts once \Gf^\ < (2561ogn/e)^ in which case we naively compute 
d{Gin,Gin,Gf^) using APSP in time 0{\Gf^\'^). Recall that even at this final point, the distances 
between marked vertices still obey the invariant. 

2.4 Running time 

We now examine the total running time of our algorithm. Let n denote the number of vertices in our 
original graph G and let V{G) denote the vertex set of the graph G in the current invocation of the 
recursive algorithm. The current invocation approximates d{Gin,Gout,Gt) as shown in subsection 
2.2 in time 0{\V {G)\/ e + 2^^^/^^). It then goes on to construct the subgraphs and G~^^^ as 
shown in subsection 2.3, where we have that after contraction using dense portals, 1 1^(6*^)1 = 
a\V{G)\+0{\og^n/e^) and |y(G+j)| = I3\V{G)\ + 0{\og^ n/ e^), where a,/3 < 2/3 and a + /3 < 1. 
The time to construct |1^(G^)| and |1^(G^„^)| is dominated by the time required to compute SSSP 
for each dense portal, which requires 0(|y(G)| • logn/e). We then continue recursively to Gf^ and 
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to Hence, if T{\V{G)\) denotes the running time for G then we get that 

T{\V{G)\) = 0{\V{G))\ ■ logn/e + 2^(1/^)) 
+ T{a\V{G)\+0{log^n/e^)) 
+ T{/3\V{G)\+Oilog'n/e^)). 

In the recursion's halting condition, once we get to components of size |y(G)| = (256 logn/e)^, we 
naively run APSP. This takes 0(|y(G)p) time for each such component, and there are 0{n/\V{G)\) 
such components, so the total time is 0(n • |y(G)|) = 0{nlog^ n/e^). It follows that 

T(n) = 0(nlog^n/e^ + n-2^(^/^)). 

3 Concluding Remarks 

We presented the first (1 + e)-factor approximation algorithm for the diameter of an undirected 
planar graph with non-negative edge lengths. Moreover, it is the first algorithm that provides a 
nontrivial (i.e. less than 2-factor) approximation in near-linear time. 

It might still be possible to slightly improve the running time of our algorithm by removing 
a logarithmic factor, or by replacing the exponential dependency on e with a polynomial one. In 
addition, the technique of Abraham and Gavoille [1] which generalizes shortest-path separators to 
the class of H-minor free graphs may also turn out to be useful. 
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